Evidence for entanglement at high temperatures in an engineered molecular magnet 
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The molecular compound [Fe2(/i2-oxo)(C3H4N2)6(C204)2] was designed and synthesized for the 
first time and its structure determined using single crystal x-ray diffraction. The magnetic suscep- 
tibility of this compound was measured from 2 to 300 K. Analysis of the susceptibility data using 
protocols developed for other spin singlet ground state systems indicates the existence of entan- 
O-l gled quantum states up to temperatures as high as 732 K, the highest entanglement temperature 

reported to date. The large gap between the ground state and the first excited state (282 K) sug- 
gests that the spin system may be somewhat immune to decohering mechanisms. Our measurements 
strongly suggest that molecular magnets are promising candidate platforms for quantum information 
processing. 
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I. INTRODUCTION 



Quantum Mechanics present unusual features, without classical analogues, such as nonlocality. Some years after 
the statement of Einstein, Podolski and Rosen on the incompleteness of Quantum Mechanics 1 -, and the possibility 
of a "spooky action at distance", Bell proposed an inequality able (in principle) to quantify all the nonclassicality 
of the correlations 2 ^, which was subsequently experimentally tested 4 -. Later, Werner showed that these quantum 
correlations cannot always be detected by this inequality 5 . Seeking a complete understanding and quantification of 
such quantum correlations, named quantum entanglement^, a whole research area was created^. 

Entanglement plays an important role in Quantum Information Science because it can be used as resource for 
many quantum information processing protocols 9 . Therefore, determining its existence in a given system is a very 
important task. From a theoretical point of view, an entangled state can be defined as a state whose density matrix 
cannot be written as a convex sum of product states^. Following this definition, some measures of entanglement could 
be constructed, i.e. entanglement of formation^ or negativity-ii, and several necessary and sufficient conditions exist 
that allow the quantification of these quantum correlations 1 ^. 

Nevertheless, from the experimental point of view, quantification of entanglement is a rather complicated task, due to 
the limited information one can obtain about the system. In some cases, experimentalists need to determine the density 
matrix of the system and then calculate the quantity of interest, which makes the process very demanding 13-15 . An 
alternative way is not to quantify the degree of entanglement in the system but only verify its existence, for example, 
using an entanglement witness^ 6 - - — (EW) . Such witnesses are observables with positive expectation value for all 
separable states and thus a negative expectation value indicates the presence of entanglement - — . This approach is 
extremely useful since one can relate these observables to experimental ones, directly detecting the contribution of 
entangled states to the measuremen t 25 ! 26 . 

In many systems, entanglement is a very fragile correlation which can vanish due to the coupling between the 
system and the reservoir, even if this coupling is weakSI. Due to this fact, it was thought that entanglement could 
only exist in low-temperature few particle systems. However, recently, it has been shown that entanglement can also 
exist in systems containing a large number of particles at finite temperatures 2 ^ - — . In this sense, the EW has helped 
to demonstrate, through the measurement of some thermodynamic observables, the existence of entangled states in 
thermal systems even at high temperature^ - — . 

Therefore, the study of entanglement in solid state physics is of great relevance to the area of quantum information 
and quantum computation, since many proposals of quantum chips are solid state based 2 ^. Furthermore, the demon- 
stration that entanglement can change the thermodynamic properties of solids, such as magnetic susceptibilit y 25 ' 42 , 
shows that entanglement can be related to significant macroscopic effects. Hence, this subject establishes an interest- 
ing connection between quantum information theory and condensed matter physics, because magnetization 4 ^, heat 
capacity^ 4 -, and internal energy 4 ^ can also be used to reveal spin entanglement among constituents of a solid. 

Quntum entanglement at elevated temperatures has been studied in several physical systems and the recent work 
by Vedral* 6 - summarizes this scenario. In that work it is recognized that some molecular magnets remain entangled 
at surprisingly high temperatures. Besides presenting robust entangled states, such materials can be engineered to 
enhance their quantum features and, consequently, to better promote the quantum information task. As an example, 
it was recently shown through magnetic susceptibility measurements, that certain families of molecular magnets are 
entangled up to 630 K, and have the maximum entanglement temperature as high as 100 K22r— . 

In this paper, we report the synthesis, crystal structure and magnetic susceptibility of a new molecular magnet: 
[Fe2(/i2-oxo)(C3H 4 N 2 )6(C204)2]. Analysis of susceptibility data suggests the existence of entangled quantum states 
up to temperatures as high as 732 K, the highest temperature reported to date. 



Recently it was showed that the magnetic susceptibility can be used as an entanglement witness, with the advantage 
of not depending on specific hamiltonians and the system can be composed of N spins S2£ Form the case of a system 
where the total hamiltonian is given by H = T-Lq+TLi, where "Ho can be any spin hamiltonian and TLi — g^BB Y)^-i S\, 
being B the weak applied magnetic field along z— direction, the only constrain is that [Ho, Hi] = 0. Thus, considering 
this constrain, the magnetic susceptibility can be written in a general form as: 



II. THERMAL ENTANGLEMENT WITNESSING 
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where g is the Lande factor, fi B is the Bohr magneton, and A 2 M Q is the variance of the magnetization in the 
a— direction (a = x, y, z). 

It shall be noted that, for a arbitrary state of a spin— 5 particle, one has^: 

{S 2 X ) + (S 2 ) + {S 2 ) = 5(5+1) (2) 

and 

(S x f + (S y ) 2 + (S z ) 2 < 5 2 (3) 

The first relation cames from the fact that the eigenvalue of the operator 5 2 is 5(5 + 1), and the second that the 
projection of the spin in any direction cannot be higher the 5. Taking the difference between these two equations it 
is found that: 

A 2 5 = A 2 S X + A 2 S y + A 2 S Z > 5 (4) 

Meanwhile, it is well known that for any separable thermal state of N— spins, one shall write that the density 
operator of the system is a convex sum of product states: 

p=Y.p-pi ) ®p2 ) ®---®p [ s ) (5) 

n 

And if it is substituted into: 

X = X* + Xv + X* = (A 2 M X + A 2 M y + A 2 M Z ) (6) 



Together with Eq.(j4]) one find that: 



k B T 



X = E [(A 2 5,)! n) + (A 2 5,); n) + (A 2 S Z ) ( : ) 



1=1 

\2 



n 

>^f^ (7) 
k B T 



In the case of isotropic systems, i.e. Xx — Xy — Xz, it is found that: 

X> ( -^f^ (8) 
3 k B i 

and one can say that the susceptibility is an entanglement witness because separable states cannot violate the relation 
showed in Eq.©. Thus, if one define: 

EW = iksT * 1 (9) 
(g^ B ) 2 NS [ > 

then for EW < there must be entangled states. 

This entanglement witness is independent of the hamiltonian of the system. It is of great importance because it 
gives the connection between the magnetic susceptibility measured in the laboratory and the existance of entangled 
states in the spin states of some magnetic material. Also, from Eq.© it is possible to determine a temperature of 
entanglement T e , i.e., the one for which the EW — ¥ and consequently that indicates the entanglement-separability 
frontier of the system^. 

To obtain T e , let us consider a Hamiltonian of a spin— 5 dimer (N = 2) composed of a Heisenberg and Zeeman 
terms: 

U = -J8 1 -S 2 -g»BB(S 1 , z + S 2iZ ) (10) 

where J is the exchange interaction, and B is an external magnetic field. Such Hamiltonian has been described 
quantitatively in the literature by the Heisenberg-Dirac-Van-Vleck Hamiltonian—^. 
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Label State \S,m s ) Dege. (B = 0) Energy (B = 0) 



E 5 \5,ms) 


11 


-15J 


E 4 |4, m s ) 


9 


-10J 


E 3 |3, ms) 


7 


-6J 


E 2 \2,m s ) 


5 


-3J 


Ei \l,ms) 


3 


—J 


Eo [0, ms) 


1 






TABLE I. States on the basis |S, ms), degeneracy for zero magnetic field and energy eigenvalues for a S — Si — S% = 5/2 
dimer. Note S ranges from \Si — S2 1 up to \Si -j- S%\. 



The magnetic susceptibility for a S = Si = S2 = 5/2 dimer and B — > ' is given by 



where 



Z = l + 3e-^ E ' + 5e"^ + 7e~^ + 9e~^ + lle^ 5 

and j3 = l/fc^T. The eigenvalues of energy Ei are in Table HI As said before, these expression describes the suscep- 
tibility of spin dimers with S — Si ~ S2 = 5/2. However, for any S = Si = S2 value, ranging from S = 1/2 up to 
S — 5/2, minor changes should be done in Eq. fTTj) . For 5 = 5/2—1/2 = 2, one should suppress the last terms of the 
numerator and of the denominator. For S = 2 — 1/2 = 3/2, one should suppress again the other last terms, and so 
on, until 5=1/2. 

Substituting Eq. fTTj) in Eq.© and calculating the boundary between separable and entangled states, which means 
EW = 0, we obtain a general form for T e = T e (J, S): 

T e = -^(2S+l)J (12) 

It should be noted that a physical temperature of entanglement (T e > 0) results only if the exchange parameter J is 
negative. 

From Eq. (fT2")) , to obtain the most robust entangled states by increasing the entanglement-separability frontier, 
one should search for the highest possible temperature of entanglement, maximizing both, the spin S value and the 
antiferromagnetic exchange parameter. To maximize S, a high spin (HS) c? 5 ion is ideal. For transitions metals the 
possibilities are either Fe(III) or Mn(II). On the other hand, to obtain a negative exchange parameter J, following 
the Goodenough-Kanamori rules, the antiferromagnetic arrangement is obtained when the angle of the d—p-d orbitals 
are close to 180°. A search in CCSD database^ revealed that there are 95 di-iron structures with Fe(III) and an 
angle Fe-O-Fe of 180°. In contrast, only 6 di-manganese structures exhibit Mn(II) and an angle Mn-O-Mn of 180°. 
Keeping these observations in mind, we exploited the syntheses of novel iron complexes using imidazole as an organic 
ligand. The present study describes a novel dinuclear Fe(III) complex which meets all of the requirements for a high 
entanglement temperature. 



III. CRYSTAL STRUCTURE 



The molecular structure of di- nuclear iron complex is presented in Fig 1(a) The crystal structure consists of a 



di-nuclear /i-oxo Fe(III) complex where the bridging oxo atom is at the crystallographic inversion center and the 
metal-to- metal distance is of 3.594(9) A with a Fe-O-Fe angle of 180°. The Fe(III) center is six-coordinated in a 
distorted octahedra coordination in which the equatorial plane is composed of three nitrogen donors from imidazole 
ligands (2.095(2), 2.147(2) and 2.168(2) A), and one oxygen atom of oxalate ligand (2.042(2) A). The apical positions 
are occupied by the /i-oxo bridging atom (1.797(5) A) and by the other oxalate oxygen (2.152(2) A), as shown in 
Fig 1(b) These Fe-0 and Fe-N distances are in the range found for other /i-oxo di-iron(III) complexes-52,. Each 
di-nuclear iron unit is linked along the b direction through N-H- • • O hydrogen bonds, between nitrogen atoms from 
imidazole ligands and the oxalate oxygen atoms, with N- ■ ■ O distances between 2.858(3) and 3.081(10) A and N- 
H- ■ • O corresponding angles ranging from 122° to 169°, forming one-dimensional (ID) network, effectively isolating 
the iron dimers (Figjl(c)'). 
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(c)Crystal packing diagram along the b axis. The hydrogen 
bonds are drawn as gray dashed lines. 



FIG. 1. Crystal structure details of [Fea^a-oxoXCaH^Na^CaO^a]. 



IV. MATERIAL DESIGN 

As can be seen from Eq.©, EW is a function of the magnetic susceptibility times temperature, \I ', of the 
material and, therefore, we measured this quantity to probe EW, as presented in Figj2](open circles - left axis). The 
entanglement witness is then easily obtained and is presented also in Fig 12 (open circles - right axis). Up to 300 K, 
EW remains negative, indicating that the system occupies entangled states up to room temperature. 

As discussed before, to maximize T e we need to maximize both, J and S (see Eq . (fT2"j) ) . To verify the HS state 
of iron in this compound, we can use the magnetic susceptibility. The Curie constant for a free 5 = 5/2 ion is 
C(S = 5/2) = 7.84 x 10" 4 ^sK/Oe, while for a S = 1/2 free ion, it reads as C(S = 1/2) = 0.67 x 10~ 4 ^sK/Oe. This 
last value is represented by an arrow in Fig |5] multiplied by two, since this figure represents the magnetic susceptibility 
for the dimer. Using this analysis it can be seen that the Fe(III) is in the HS - S — 5/2 configuration. 

Before fitting Eq. (|lip to the experimental data, it is noticeable that the quantity yT does not approach to zero 
when the temperature goes to zero (see Fig[2]). this background contribution arises from a quite small amount p of 
unreacted free Fe(III) ions. This amount does not affect the discussion of the present paper, since our aim is hundreds 
of Kelvin above room temperature. Thus, the final expression to be fitted to the experimental data is: 

2C(S = 5/2) 

X=(l-p)Xd + P T (13) 
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FIG. 2. Left axis: Experimental (open circles), and theoretical (solid lines) magnetic susceptibility times temperature for the 
S — Si — S2 — 5/2 Iron dimer. The thicker red line represents the fitting of Eq. (|13|l to the data; while the thinner black line is 
an extrapolation, using the optimized parameters, up to higher temperatures. Right axis: Entanglement witness as a function 
of temperature. Open circles are from the experimental data, while the solid line is an extrapolation up to higher temperatures 
using the theoretical model (Eqs.© and (|13|> ) and the optimized parameter (see text). Note that this material has entangled 
states up to 732 K. 



where C(S = 5/2) = 7.84 x 1CT 4 /Lt B K/Oe is the Curie constant of S = 5/2 free ion. The red thicker line in Figgis 
the fitting of the above model to the experimental data. The optimized parameters are: g = 1.8, J = —282 K and 
p = 1.0%. Note that the first term of the above equation has a factor of 2, because the susceptibility is expressed for 
a dimer. Based on the results above, it is possible to extrapolate the susceptibility curve up to higher temperatures, 
as can be seen in Figj2] (thiner black line). This quantity in Fig{2] tends to the Curie constant of a free S = 5/2 ion. 

Thus, the EW emerges and we can verify in Fig that this Iron dimer has entangled states that can survive up to 
732 K, hundreds of Kelvin above room temperature, and well above the expected thermal stability of this compound. 
This result is in accordance with the empirical prediction of Eq. (IT2l) which gives T e — 761 K ( see Ref.[4l|). Therefore 
this compound shows the highest entanglement temperature reported in the literature, for solid state systems^ 6 -. 



V. CONCLUSION 



The robust entangled states exhibited by molecular magnets make them among the most promising materials for 
quantum information processing applications. 46 . Based on our previous work, to achieve the highest entanglement 
temperature promising compounds should contain antiferromagnetically coupled magnetic ions with the highest pos- 
sible spin states 41 . A systematic study of molecular iron compounds was made, which resulted in the preparation of 
a new di-iron complex with formula [Fe2(/i2-oxo)(C3H4N2)6(C204)2]. Analysis of magnetic susceptibility data from 
this compound indicated quantum entanglement (or the entanglement-separability frontier) up to temperatures as 
high T e = 732K. This is the highest entanglement temperature reported to date^ -. 



Appendix A: Synthesis 



All reagents and chemicals were purchased from commercial sources and used without further purification. Iron 
oxalate dihydrate (99%) and imidazole (99%) were obtained from Aldrich and Panreac Chemical, respectively. A 
mixture of iron oxalate dihydrate (1.39 g; 7.90 mmol), imidazole (0.52 g; 7.90 mmol) and ethanol (70ml), were 
refluxed for three hours. The resultant brow-orange solution was allowed to cool down overnight. A light brown 
coloured crystalline product was obtained after 3 weeks. 



Appendix B: Crystallography 



The newly synthesized Fe(III) complex has been structurally characterized by single crystal X-ray crystallography. 
The complex with formula [Fe2(/J.2-oxo)(C3H4N2)e(C204)2] crystallizes in the triclinic crystal system with the space 
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Bond lenghts (A) 



Fe-O(l) 
Fe-N(8) 


1.797(5) 
2.147(2) 


Fe-0(2) 2.152(2) 
Fe-N(13) 2.095(2) 


Fe-0(3) 
Fe-N(18) 


2.042(23) 
2.168(2) 


Bond angles (°) 


0(l)-Fe-0(2) 

0(1)-Fe-N(8) 

0(2)-Fe-N(8) 

0(3)-Fe-N(8) 

N(8)-Fe-N(13) 

Fe-0(l)-Fe A 


175.26(5) 

96.56(6) 

81.90(8) 

89.09(8) 

92.29(9) 

180 


0(l)-Fe-0(3) 97.41(5) 
0(1)-Fe-N(13) 97.47(6) 
0(2)-Fe-N(13) 87.09(7) 
0(3)-Fe-N(13) 164.80(7) 
N(8)-Fe-N(18) 164.59(8) 


0(2)-Fe-0(3) 78.12(7) 
0(1)-Fe-N(18) 98.79(6) 
0(2)-Fe-N(18) 82.69(8) 
0(3)-Fe-N(18) 87.55(8) 
N(13)-Fe-N(18) 87.12(9) 



(A) 1 — x, 2 — y, 1 — z. 



TABLE II. Selected bond lengths (A) and angles (°) for [FeaQte-oxoXCa^NaMCaO^a]. 



D-H- • ■ A 


H---A (A) 


D---A (A) 


D-H- • ■ A (A) 


N(ll)-H(ll)- 


•0(7) [-l+x,y 


A 




2.05 


2.922(4) 


169 


N(16)-H(16)- 


•0(5) [x,l + y,z 






2.45 


3.013(3) 


122 


N(16)-H(16)- 


•0(7) [x,l+y,z 






2.01 


2.858(3) 


161 


N(21)-H(21)- 


•0(3) [2-a;,2- 


yA- 


-*] 


2.48 


3.081(10) 


126 


N(21)-H(21)-- 


■0(5) [2-31,2- 




-*] 


2.17 


3.028(12) 


164 



TABLE III. Hydrogen bond dimensions for [Fea^a-oxoXCs^Na^CaO^a]. 



group PI where a = 8.4789(11), b = 9.6481(12), c = 9.769(2) A, a = 104.350(7), /3 = 91.673(7), 7 = 115.215(5) 
°,V = 692.1(2) A 3 , Z = 1. The X-ray data was collected on a CCD Bruker APEX II at 150(2) K using graphite 
monochromatized Mo-Ka radiation (A = 0.71073 A). The crystal was positioned at 35 mm from the detector and the 
frames were measured using a counting time of 60 s. At the total, 7211 reflections were collected and subsequently 
merged to 3383 unique reflections with a i?mt of 0.0374. Data reduction including a multi-scan absorption correction 
were carried out using the SAINT-NT from Bruker AXS. The structure was solved using SHELXS-97^ 3 . and refined 
using full-matrix least squares in SHELXL-97 5 — . The carbon and nitrogen bonded hydrogen atoms were included at 
calculated positions. Anisotropic thermal parameters were used for all non-hydrogen atoms while the hydrogen atoms 
were refined with isotropic parameters equivalent to 1.2 times those of the atom to which they were attached. Two 
atoms N(21) and C(22) from one imidazole group were found to be disordered over two positions, and these were 
introduced in the structure refinement with adjustable occupancies of x and 1 — x, x being equal to 0.568. The final 
refinement of 225 parameters converged to final R and R w indices R\ — 0.0411 and WR2 = 0.0856 for 2495 reflections 
with / > 2 (I) and i?i = 0.0683 and V0R2 = 0.0954 for all hkl data. The molecular diagrams presented are drawn with 
graphical package software Diamond. The relevant bond distances and angles are reported in Table [TT] The detailed 
hydrogen bond distances and angles are reported in Table IIIII 

The crystallographic data for [Fe2(/J.2-oxo)(C3H 4 N 2 )6(C204)2] complex has been deposited in Cambridge Crys- 
tallographic Data Centre (CCDC 879631). The data can be obtained free of charge at www.ccdc.cam.ac.uk or from 
Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, UK; fax +44 (0) 1223 336033, e-mail: 
deposit@ccdc.cam.ac.uk. 
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